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ABSTRACT 

Context. We present here the results of a 180 ks Ciiandi'a-LETGS observation as part of a large multi- wavelength campaign on Mrk 

509. 

Aims. We study the warm absorber in Mrk 509 and use the data from a simultaneous HST-COS observation in order to assess whether 

the gas responsible for the UV and X-ray absorption are the same. 

Methods. We analyzed the LETGS X-ray spectrum of Mrk 509 using the SPEX fitting package. 

Results. We detect several absorption features originating in the ionized absorber of the source, along with resolved emission lines 

and radiative recombination continua. The absorption features belong to ions with, at least, three distinct ionization degrees. The 

lowest ionized component is slightly redshifted (Av = 73 km s"') and is not in pressure equilibrium with the others, and therefore it is 

not likely part of the outflow, possibly belonging to the interstellar medium of the host galaxy. The other components are outflowing at 

velocities of -196 and -455 km s"', respectively. The source was observed simultaneously with HST-COS, finding 13 UV kinematic 

components. At least three of them can be kinematically associated with the observed X-ray components. Based on the HST-COS 

results and a previous FUSE observation, we find evidence that the UV absorbing gas might be co-located with the X-ray absorbing 

gas and belong to the same structure. 

Key words. Galaxies: individual: Mrk 509 - Galaxies: Seyfert - quasars: absorption lines ~ X-rays: galaxies 

1 . Introduction Interestingly, sources that present absorption in the UV also 

show absorption in X-rays (Crenshaw et al. 1999), which may 
point to a possible relation between the two phenomena. With 

. ^. , ^. , ■ / A ^xTx 1 1 • the advent of modern medium- and high-resolution spectro- 

Active galactic nuclei (AGN) are among the most luminous per- ..... ,. .° . ,, 

. ^ ° ^ , . , • .u TT • A • vc . graphs, the ionization state, kinematic properties and column 

sistentextragalactic sources known in the Universe. A significant j.- ^, , ■ \c ■ <t uij 

„ . J- , • , , • , • • • ■ J ^r densities or the absorption features in X-rays could be deter- 

rraction or their bolometric luminosity is emitted at X-ray wave- • j • i i • i 1 1 , , , 

, , . ■ ■ ■ , ■ ■ r , mined with higher accuracy. However, even though the average 

lengths, and originates in the innermost regions or these sources, ^ , .° j • i -i^ • n • 

,° ■ 1.1 1 u 1 J ,^ •. .• 1 .• outflow velocity measured in the X-rays IS generally consistent 

where a supermassive black hole resides. Gravitational accretion ... i • , •-• j • ■ rrir i i 

„ ^^ ^ .1 • 1 1 1 1 1 ■ 1 1- J 1 .1 1 • with the velocity shirts measured in the UV, the spectral reso- 

or matter onto this black hole IS believed to be the mechanism ,. .,, -^ . . „. , , j-ji 

, ^ , , . ^ J .1 • ^ J- ■ z: ij 1 lution in X-rays is insumcient to resolve the many individual 

that powers these objects, and the intense radiation held thus , . ^. -^ ^ ^- ^i ttat 

^ J ■ u. J ■ .a c .u 1 Tjj kinematic components present in the UV. 

generated might drive gas outflows from the nucleus. Indeed, ^ ^ 

more than 50% of the Seyfert 1 galaxies presents evidence of a Some authors discussed the possibility that at least a fraction 

photoionized gas, the so-called warm absorber (WA hereafter), of the UV absorption is produced in the same gas responsible for 

in flieir soft X-ray spectra (e.g. Reynolds [1997] George et al. the X-ray absorption (Mathur et al l 1 9941 [19951 1 1 999t Crenshaw 

[1998 Crenshaw et aI. l2003aD . These WA are usually blueshifted & Ki-aemer 1999 Kiiss et al. [20001 Kraemer et al. "2002), but 



with respect to the systemic velocity of the source with velocities the results were inconclusive and the connection between the 
of the order of hundreds of km s ' (Kaastra et al. 120001 Kaspi et two absorbing components seems to be very complex. The most 
al. 1200 1 1 Kaastra et al. 120021 ). common source of uncertainty came from the lack of correspon- 
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dence between the column densities of the species detected in 
both UV and X-rays. However, more recent high-resolution ob- 
servations indicate that these values can be reconciled if one as- 
sumes a velocity-dependent partial covering in the UV (e.g. Arav 
etal. 120021 120^ . 

Simultaneous UV/X-ray observations are therefore the key 
to draw a global and unbiased picture of the WA in AGN (see e.g. 
Costantini 2010 for a review). High-resolution UV observations 
covering a wide range of ionization states are essential to deter- 
mine whether some of the UV absorbers are indeed associated 
with the X-ray warm absorber Low-ionization absorbers, such 
as Mg II and C ii, are only detected in the UV domain, while high- 
ionization ions like O viii or Mg xi can only be detected in X- 
rays. However, there is a range of ionization states that may pro- 
duce ions detected in both bands, namely N v and, particularly, 
O VI. Therefore, several simultaneous or quasi-simultaneous X- 
ray and UV spectroscopic observations have been carried out 
on many sources known to host a X-ray WA and UV absorbers, 
e.g. NGC 3783 (Kaspie t al. l2002l Gabel e tal.l2003T l. NGC 5548 
(Crenshaw et al. '2003b' Steenbrugge et al.'2005^, Mrk 279 (Arav 
et al. 2005, Gabel et al. 2005 Costantini et al. 2007). 

Mrk 509 is a Seyfert 1/QSO hybrid and is one of the 
best studied local AGN because of its high luminosity (L(l- 
1000 Ryd)^ 3.2 x 10"^^ erg s"') and proximity. It is particu- 
larly suitable for an extensive multi-wavelength campaign be- 
cause of its brightness and confirmed presence of an intrinsic 
WA (Pounds et al. 120011 Yaqoob et al. l2003l Smith et al. l2007l 
Detmers et al 



in X-rays, as well as a variety of UV absorp- 
tion li nes (C renshaw et al. 1995 1999; Savage et al. ll997l Kriss 
et al. 120001 Kraemer et al. 2003 ) and slow variability, which 
makes it particularly suited to reverberation studies. 

This paper is devoted to the analysis of the Qiandra-LETGS 
observation of Mrk 509, which is part of an unprecedented multi- 
wavelength campaign of observations carried out at the end of 
2009, aiming at determining the location of the warm absorber 
and the physics of its ionized outflow. The campaign involved 
five satellites (XMM-Newton, Chandra, INTEGRAL, HST and 
Swift), and two ground-based facilities (WHT and PAIRITEL), 
thus covering more than five orders of magnitude in frequency. 
An overview of the campaign and its goals can be found in 
Kaastra et al. ( 1201 lal ). hereafter Paper I, and references therein. 

This paper is organised as follows: in Sect. |2] we briefly 
present the reduction method applied, in Sect. |3] we describe 
the different models used to fit the spectrum, while in Sect. |4] 
we discuss our results and compare them with previous X-ray 
observations as well as with the simultaneous HST-COS obser- 
vation. Finally, our conclusions are summarized in Sect. |5] The 
quoted errors refer to 68.3% confidence level (A^^ = 1 for one 
parameter of interest) unless otherwise stated. 



2. X-ray observations and data reduction 

Mrk 509 was observed on two consecutive orbits in December 
2009 by the Low Energy Transmission Grating Spectrometer 
(LETGS, Brinkman et al. '2000^ coupled to the HRC-S detec- 
tor on board of Chandra for a total exposure time of ~180 ks 
(see Table[T]). The LETGS observation was made simultaneously 
with the HST-COS spectrograph observations (Kriss et al. 1201 II) . 
The LETGS data were reduced using the standard pipeline 
(CIAO 4.2) until the level 1.5 event files were created. The rest 
of the procedure, until the creation of the final level 2 event file, 
followed the same steps as the standard pipeline regarding the 
wavelength accuracy determination and effective area genera- 
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Fig. 1. Lightcurve of the Chandra observations. The solid hori- 
zontal lines mark the HST-COS observation intervals. 

Table 1. Qiandra-LETGS observation log. 



Obs ID Date 

(dd-mm-yyyy) 



Exp. time Count rate 
(ks) (counts/s) 



11387 
11388 



10-12-2009 
12-12-2009 



131.4 
48.6 



0.94 
1.05 



tion, but was carried out using an independent procedure first 
described in Kaastra et al. ( 120021) . 

The lightcurve of both observations is shown in Fig.[T] with 
a bin size of one hour (3600 s). We overplot the time intervals 
when the HST-COS observations were taken. The lightcurve 
shows no episodes of short-time variability (i.e. within a few 
kiloseconds) when we apply a bin size of 100 s. At longer 
timescales the source flux varies smoothly throughout the whole 
observation. There seems to be a drop in flux, starting shortly 
after the beginning of the observation, followed by a recovery at 
f ~ 50 ks. The deviation from the average flux was of the or- 
der of 30% and no significant spectral variations were observed 
between the higher and lower flux states. 



3. Spectral analysis 

The combined LETGS spectrum was analysed using the SPE5U 
fitting package v2.03.00 (Kaastra et al. |1996l l. The spectrum has 
variable bin sizes, depending on the statistics, but is typically 
binned in 0.03 A bins, which results in ~30 counts per bin. 
Therefore, the fitting method adopted was x^ statistics. In what 
follows we describe the different models used to describe the 
continuum, and the absorption and emission features present in 
the spectrum of Mrk 509. The cosmological redshift of Mrk 509 
is z = 0.034397 (Huchra et al. .l993j . For this work, we adopt a 
systemic redshift of z = 0.03448, which is the heliocentric red- 
shift of Mrk 509 transferred to the Earth-centered frame of our 
Chandra data at the time of our observation. Note that the stan- 
dard CIAO software does not take this motion into account. An 
overall view of the LETGS spectrum of Mrk 509 in the range 
where the majority of relevant lines fall is shown in Fig.|2] 



' |http : //www ■ sron . nl/spex] 
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Table 2. Best-fit spline continuum parameters. 



20 25 

Wavelength (A) 

Fig. 2. Chandra-LETGS spectrum in the 9 - 38 A range. The 
continuum is described with a spline model. The most prominent 
lines are labeled. 



Wavelength (A) 


Flux" 


4.00 







4.86 


11.5 


±8.0 


5.90 


11.3 


±0.1 


7.16 


10.8 


±0.2 


8.69 


11.0 


±0.1 


10.55 


11.1 


±0.2 


12.81 


12.0 


±0.2 


15.55 


12.8 


±0.2 


18.87 


14.0 


±0.3 


22.91 


15.6 


±0.3 


27.82 


17.3 


±0.2 


33.78 


19.3 


±0.3 


42.42 








Notes. '"' Unabsorbed flux, in units of 10 "* photons cm ^ s 'A '. 

Table 3. Best-fit power-law plus modified black body continuum 
parameters. 



3.1. Continuum 

The continuum was modelled following two different ap- 
proaches. We first used a spline model, that reproduces the shape 
of the continuum without any a priori assumptions on its physi- 
cal origin. The boundaries of the spline were set between 4 and 
41 A. We note that the analysis of the spectrum is nominally car- 
ried out in the range 5 - 40 A, where the majority of the most 
relevant spectral features fall. The spline boundaries had to be 
extended beyond these limits in order to prevent spurious arti- 
facts in the proximity of the boundaries, since the model sets the 
photon spectrum to zero outside the defined limits. The grid con- 
sisted of 13 points spaced logarithmically in wavelength. This 
spacing is wide enough to successfully model the continuum 
without risking to fit any intrinsic absorption features. The best- 
fit flux at the hinge points of the grid is reported in Table |2] 

As an alternative approach, we also modelled the continuum 
using a single power-law in the same range as the spline model. 
This yielded an unacceptable fit, both visually and in terms of 
reduced;^^, due to deviations at soft X-rays. To account for these 
deviations we added a modified black body model (mbb model 
in SPEX). The final fit yielded a comparable reduced ^^ (1742 
for 1329 d.o.f.) to that of the spHne model (1721 for 1320 d.o.f.). 
The best-fit parameters are reported in Table|3] In both cases, the 
continuum was corrected for local foreground absorption (see 
SectO. 

The spline and power-law plus mbb unabsorbed models are 
compared in Fig. [3] where the rest of emission components are 
also included. Both models agree well for most part of the fit- 
ting range, with average deviations between both models of less 
than 5% in the 10 - 36 A range. The larger deviations occur 
below 10 A and above 36 A (where the maximal deviations 
are only ~7% and ~9%, respectively), most likely due to the 
logarithmic sampling of the spline model. Indeed, the power- 
law model fails to correctly reproduce the continuum at shorter 
wavelengths, which is best covered by the spline with 4 hinge 
points below 10 A. The poorer sampling of the spline at longer 
wavelengths with only one hinge point above 30 A provides a 
slightly worse description of the continuum than the power-law 
plus mbb model. However, these deviations are lower or similar 
to the statistical uncertainties in the data, which are of the order 
of 10%, and therefore the continuum models do not critically dif- 



Model 



'ra 



-fo.5 



Power law 2.06i 
Black body 



:0.03 



0.126±0.005 



2.56±0.05 
0.22±0.03 



Notes. '"' Temperature, in units of keV: 
<*' 0.5-2 keV flux, in units of 10"" 



erg cm ^ s ' . 



fer. Since the differences between both models in terms of good- 
ness of fit are better for the spline model by A;\f/Av = 21/9, this 
was the baseUne model adopted for the fit in the 5 - 40 A range. 
The extension of the analysis towards longer wavelengths, up to 
~60 A, was made using only the power-law baseline model (see 
Sect. Ill- 



3.2. Local absorption 

Before analyzing the spectral features intrinsic to Mrk 509 we 
first account for the foreground absorption of our Galaxy. The 
local interstellar medium (ISM) presents three different phases: 
a cold neutral phase, a warm slightly ionized phase, and a hot 
highly ionized phase (Pinto et al. 120101) . For the neutral phase 



we take Nh = 4.44 x 10^" cm^-^ (Murphy et al.TWS). The local 
absorption parameters were kept fixed during the spectral fits. A 
detailed analysis of the local ISM absorption in the line of sight 
of Mrk 509 will be reported in a forthcoming paper (Pinto et al., 
in preparation). 



3.3. Emission lines 

There are several emission features that can be visually identi- 
fied in the LETGS spectrum of Mrk 509. In most of the cases 
the emission lines are blended with features from the warm ab- 
sorber. They need, however, to be modelled in order to account 
for the excesses seen on both sides of the absorption line. This 
was the case of the C vi Lye line at ~34.9 A and, with less sig- 
nificance, the Oviii Lya and Nvii Lya lines at ~19.6 A and 
~25.5 A, respectively. The Ovii region shows also great com- 
plexity, as a number of emission lines are blended with Ovi 
and O VII absorption lines. Since the lines look broad, we have 
modelled the O vii triplet, as well as the lines mentioned above, 
with Gaussians. The high flux ratio of the O vii intercombination 
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20 30 

Wavelength (A) 



40 



Fig. 3. Unabsorbed continuum plus emission line models. Solid 
line: spline continuum; dot-dashed line: power-law + modified 
black body. 

Table 4. Best-fit emission line parameters. 



Line 


Wavelength" 


FWHM" 


Flux'- 


Ovii(r) 


21.602 


0.285 


<0.43 


vii (i) 


21.808 


0.288 


0.52±0.20 


Ovn(f) 


22.101 


0.292 


0.67±0.20 


N VII Lya 


24.781 


0.327 


<0.42 


VIII Lyo- 


18.973 


0.250 


<0.31 


CviLyo' 


33.736 


0.445 


0.67±0.25 



Notes. *°' Rest wavelength in A (frozen); 
*** Full width at half maximum in A (frozen); 
*''' In units of 10""* photons cm"^ s"'. 

Table 5. Best-fit RRC parameters. 



Continuum 


Ion 


EM" 


T" 


Spline 


Cvi 

Nvii 


'•■^-4.5 
-"•■^-2.1 


3 7+3** 

■-'•'-2.0 

7 1 +'■** 


PL+mbb 


Cvi 
Nvii 


6 2+^-'^ 

4.2!^:^ 


3.4-i 
7.9-;? 



Notes. *"' Emission measure n(,nio„ V, in units of 10*' cm ^\ 
'*' Temperature, in units of eV. 



line with respect to the other lines of the O vii triplet suggests 
that the gas has high density, thus pointing towards an origin in 
the BLR. This was also found in other Seyfert 1 galaxies (e.g. 
Longinotti et al. 2008 , 2009 ), and is generally different from the 
results found in Seyfert 2 galaxies, where the lines are likely 
produced in the NaiTow Line Region (NLR, e.g. Bianchi et al. 
12006 ' Guainazzi & Bianchi l2007l l. Therefore, the wavelengths 
were fixed to the rest-frame values and the widths to those of 
the optical BLR as measured for the B aimer lines with XMM- 
Newton OM (FWHM ^ 4000 km s ', Mehdipour et al. l20m ). 
The fluxes of the lines are therefore the only free parameters. The 
O VII forbidden (/) and intercombination (/), and the C vi lines 
are detected at the ~ 3cr confidence level, whereas the others are 
just detected at the Icr confidence level. The best-fit values are 
summarised in Table |4] 




20 
Wavelength (A) 



40 



Fig. 4. Model plot showing the contribution to the spectrum of 
each of the three WA components. The solid line represents com- 
ponent 1, while the dashed and dot-dashed lines represent com- 
ponents 2, and 3, respectively . For clarity, components 2 and 3 
have been shifted in flux by a factor of 1 .5 and 2, respectively. 



In addition to these broad lines, there are also hints of radia- 
tive recombination continua (RRC) present in the spectrum. We 
found a RRC of C vi, seen at ~26.2 A, at the > Icr confidence 
level, and a marginal (detected at the 1.65cr confidence level) 
RRC of N vii around 19. 1 A. Other RRC like that of O vii were 
statistically undetectable in our dataset. The features were fitted 
using an ad-hoc model available in SPEX that takes into account 
the characteristic shape of these lines. The best-fit temperature 
and emission measure of these features are reported in Table |5] 
The best-fit low temperatures found for these features indicate 
that the dominant mechanism in the gas is photoionization. The 
emission line models used in the fit are shown in Fig.[3]for both 
the spline and power-law continuum models. It can be seen that 
the choice of either continuum modelling does not significantly 
affect the detection of the emission lines. 

3.4. Warm absorber 

The majority of the absorption lines seen in the soft X-ray spec- 
trum of Mrk 509 corresponds to the presence of ionized gas at 
the redshift of the source. Most of them are blueshifted with re- 
spect to the systemic velocity of Mrk 509. 

We first modelled the WA using two slab components to 
describe the low- and high-ionization absorption lines, respec- 
tively. The slab model in SPEX calculates the transmission of a 
slab of material without any assumptions on the underlying ion- 
ization balance. The free parameters are the ionic column densi- 
ties, the nns velocity broadening, and the outflow velocity. The 
slab fit showed that the low-ionization ions, i.e. those for which 
their concentrations peak at log ^ lower than 2.1, appear to be in 
a different kinematic regime with respect to the high-ionization 
ions. The ionization parameter is defined as ^ = L/nR^, where 
L is the ionizing luminosity in the 1 - 1000 Ryd range, « is the 
density of the gas, and R is the distance from the ionizing source. 
Ions belonging to the first slab component are at the systemic 
velocity, whereas those of component 2 are outflowing at sev- 
eral hundreds of km s ' (model 1 in Table |6]l. This model is, 
however, not entirely realistic, as it does not allow for a possible 
multi-velocity structure in the outflow. Different kinematic com- 
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Table 6. Ionic column densities of the most relevant ions. 





Model 1 




Model 2 




Component 


1 


2 


1 


2 3 




a" 


108±11 


100+«3 


85+17 


61 + 15 90+36 




Vout* 


3+17 


-338!Ti, 


87+36 


-208+42 -635+172 




Ion 


logWion'- 


logM™' 


logA'ion^ 


logiV,™,'' logM„„'' 


logf" 


Cv 


16.47+0.17 




16.38+0.19 


15.95+0.43 


0.15 


Cvi 


16.66±0.19 




16.65+0.18 


15.85+0.41 


1.20 


Nv 


<15.4 




<15.3 


<15.4 


0.00 


Nvi 


15.57+0.36 




15.60+0.37 


<15.3 


0.75 


Nvn 


15.89+0.27 




15.99+0.23 


<15.5 


1.60 


Om 


16.16+0.23 




15.96+0.36 


<15.6 


-1.50 


Oiv 


16.07+0.31 




15.76+0.31 


<15.9 


-0.65 


Ov 


16.17+0.22 




<15.6 


16.45+0.22 


-0.05 


Ovi 


16.19+0.18 




15.99+0.22 


15.69+0.39 


0.40 


Ovii 


17.10+0.11 




16.76+0.14 


16.54+0.22 


1.15 


Ovm 


17.43+0.13 




16.76+0.31 


17.14+0.18 


1.90 


Ne VIII 


<16.0 




<16.0 


<15.7 


1.20 


Neix 


16.49+0.29 




16.54+0.31 


<15.7 


1.75 


Nex 




16.64+0.27 


<15.8 


16.37+0.35 16.16+0.40 


2.40 


Mgxi 




16.77+0.42 


<16.1 


<16.2 16.90+0.33 


2.25 


Six 


15.62+0.33 




15.40+0.34 


15.39+0.38 


1.90 


Si XI 




<15.2 


<15.2 


<15.1 


2.10 


Svm 


<15.2 




<15.1 


<15.3 


0.80 


Six 


15.10+0.30 




<15.2 


15.17+0.40 


1.35 


Sx 


<15.0 




14.75+0.30 


<14.7 


1.70 


Sxi 


<15.3 




<15.5 


<14.8 


1.95 


Sxn 




<15.6 


<15.3 


<15.7 


2.20 


Sxm 




<14.9 


<15.4 


<14.8 


2.40 


Sxiv 




<15.7 


<15.8 


<15.5 


2.60 


Arix 


<15.1 




<15.0 


<15.3 


0.90 


Arx 


<15.4 




<15.3 


<15.2 


1.35 


Arxi 


<14.4 




<14.7 


<14.8 


1.70 


Arxn 


<14.5 




<14.5 


<14.3 


2.05 


Caxiv 




15.53+0.26 


<15.2 


15.57+0.24 


2.35 


Fei 


<16.0 




<15.1 


<15.4 


-8.50 


Fen 


<16.1 




<15.6 


<15.6 


-8.50 


Fern 


<16.1 




15.69+0.43 


<15.5 


-2.10 


Feiv 


<16.0 




<15.5 


<15.8 


-1.40 


Fev 


<15.3 




<15.0 


<15.0 


-0.95 


Fevi 


<15.2 




<15.2 


<15.0 


-0.55 


Fevn 


<15.3 




<15.0 


<15.0 


0.05 


Fevm 


<15.6 




<15.4 


<15.0 


0.85 


Feix 


15.44+0.27 




15.41+0.26 


<15.0 


1.40 


Fex 


<15.0 




<14.7 


<14.9 


1.65 


Fexi 


<15.5 




<15.3 


<15.1 


1.85 


Fexn 


<15.2 




<15.2 


<15.4 


2.00 


Fe xm 


<15.3 




<15.2 


<15.2 


2.05 


Fexiv 




<15.5 


<15.5 


<15.1 


2.10 


Fexv 




<15.3 


<15.0 


<15.0 


2.15 


Fexvi 




15.64+0.22 


<15.0 


15.67+0.23 


2.20 


Fe xvn 




<15.3 


<14.7 


<15.0 


2.30 


Fe xvni 




15.84+0.17 


<15.0 


15.76+0.19 


2.50 


Fexix 




<15.7 


<15.1 


<15.6 


2.80 


Fexx 




15.68+0.31 


<15.1 


15.56+0.43 


3.00 


Fexxi 




15.80+0.32 


<15.6 


<15.5 15.81+0.32 


3.20 


Fexxn 




<15.9 


15.92+0.50 


<15.3 


3.30 


Fe xxiii 




<16.2 


<16.2 


<15.8 


3.40 


Fe XXIV 




<16.5 


<16.5 


<16.3 


3.50 



Notes. *"' rms velocity broadening, in units of km s"' ; 

'*' Outflow velocity, in units of km s"'; 

''^' Ionic column density, in units of cm"^; 

'''' Ionization parameter at which the concentration of the ions peaks, in units of erg s"' 



ponents may, in principle, contribute to the total column den- 
sity of a given ion. Therefore, the model was extended and all 
the ions were fitted with both velocity components (model 2 



in Table |6]l. The results show that indeed Cv and Cvi, as well 
as O VI, O VII and O vm ions are significantly detected in both 
velocity regimes. Moreover, a third component is required in 
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order to correctly fit Mgxi and Fexxi which seem to be out- 
flowing at even larger velocities. Hence, the final fit shows three 
distinct kinematic components: one slightly redshifted and two 
blueshifed with respect to the systemic velocity of the source 
by ~200 and ~600 km s"'. The best-fit values for models 1 and 
2 are reported in Table |6] along with the ionization parameter 
^ for which the ion fraction peaks for that ion. The values of 
^ are derived from an ionization balance model and are shown 
for illustrative purposes, as the slab fit is not model-dependent. 
These different velocity components of the WA will be used as a 
starting point for a more physical analysis. 

We also use an alternative model to describe the WA in Mrk 
509. The xabs model in SPEX calculates the transmission of 
a slab of material where all ionic column densities are linked 
through a photoionization balance model, which we calculated 
using CLOUDY (Ferland et al. 1998) version COS. The spec- 
tral energy distribution (SED) provided to CLOUDY as an input 
was built using the different observations involved in this cam- 
paign (Chandra-LETGS, HST-COS, INTEGRAL, SWIFT), plus 
archival FUSE data. Details on how the SED was built and plots 
are shown in Paper I. The free parameters in the xabs fit are 
the ionization parameter ^, the hydrogen column density Nh, the 
line velocity width cr, and the outflow velocity Vo„i- The abun- 
dances were set to the proto-Solar values of Lodders & Palme 

mm . 

A single xabs component yields a^^ of 1490 for 1310 d.o.f. 
but it is insufficient to describe the different velocity compo- 
nents and the range of ionization states found in the slab fit. 
A second xabs component was added, further improving the fit 
in A;\f^/Av = 46/4. These two components with log^ ~ 1 and 
log^ ~ 2, respectively, account for the lowly and the major- 
ity of the highly ionized ions. However, a few high-ionization 
lines located at shorter wavelengths (namely Neix, Fexx, and 
Fe xxi), although weaker, remain unfitted. We therefore added a 
third xabs to our model to fit the fast (vout 460 km s ' ) high- 
ionization (log^ ~ 3) component in the WA. This component is 
detected slightly above the 2cr level (the fit improves a further 
A;t'^/Av - 18/4), giving also a visually acceptable fit of these 
ions. Adding more xabs components no longer improves the fit, 
thus we recover three kinematic components in the WA in agree- 
ment with the results of the slab fit. A plot with the contribution 
to the spectrum by each xabs component is shown in Fig.lD 

In summary, we observe in Mrk 509 three WA components 
with distinct ionization degrees (with log^ =1.06, 2.26, and 
3.15, respectively). All three WA phases are somewhat shal- 
low (with A^H a few times 10^'^ cm"^), with the first component 
slightly redshifted with respect to the systemic velocity of Mrk 
509. The second and third components are outflowing at veloci- 
ties of ~ 200 and ~ 460 km s"', respectively. The best-fit results 
are shown in Table |7] The detailed LETGS spectrum along with 
the best-fit xabs model is shown in Fig.|5] 

This WA fit was performed with a spline continuum. Since 
the continuum was also modelled following a different approach 
(power-law plus modified black body) that yielded a similar 
goodness of fit (as explained in Sect l3.1b . we repeated the WA 
analysis using this continuum model for comparison. In this fit, 
the nns velocity broadening cr was kept fixed to those values ob- 
tained in the fit using a spline as a continuum. As can be seen in 
Table|2l the WA parameters obtained are, as expected, fully con- 
sistent with those assuming a spline model for the continuum. 

Furthermore, we extended the analysis range towards longer 
wavelengths, trying to take advantange of the high effective area 
of LETGS up to ~170 A. Unfortunately, the high background 



Table 7. Warm absorber best-fit parameters. 



Continuum Comp. 



logr 



Vout 



Spline 



1 


1.06±0.13 


1.9±0.5 


S6-» 


73+41 
'■-'-59 


2 


2.26±0.07 


5.4±1.6 


49!?^ 


-mt^l 


3 


3.15±0.15 


5 9+^' 

■'■^-2.6 


75-f 


-455+'"^ 



PL+mbb 1 1.06±0.14 1.7+0.4 86 (/) 73±38 

2 2.26+0.07 5.1 + 1.3 49 (/) -202+g 

3 3.35+0.27 8.1X^-5* 75(/) -445*\ll 

Notes. *°' Ionization parameter, in units of erg s"' cm; 
**' Column density, in units of 10^" cm"^; 

*^' nns velocity broadening, in units of km s"', (/) denotes a fixed pa- 
rameter during the fit; 
'''' Outflow velocity, in units of km s"'. 



made only possible to carry out an analysis up to ~60 A. The 
most relevant ions in this part of the spectrum correspond to 
C v. Si XI, Si X, S IX, and S x. Since the bulk of the WA lines 
falls below 40 A, the WA parameters described with xabs are 
tightly constrained. Therefore, an extended fit including the 
longer wavelengths provided results fully consistent with those 
reported in tableQ A slab fit for these ions provided upper lim- 
its for the majority of them with the exception of Si x, S ix, and 
S X (see Table |6]l, albeit with large error bars. 



4. Discussion 

4. 1 . Comparison with previous X-ray observations 

Prior to this LETGS observation, Mrk 509 has been analyzed 
using various high-resolution X-ray spectra. It was observed by 
XMM-Newton in October 2000 for -30 ks. Pounds et al. (|200TT) 
performed an analysis of the RGS spectrum that showed hints of 
intrinsic absorption, mainly Ne ix features possibly blended with 
Fe XIX, as well as marginal emission features in the O vii com- 
plex. An additional ~30 ks XMM-Newton observation was made 
in April 2001. Smifli et al. ( 120071) combined both RGS datasets 
in order to increase the signal-to-noise ratio. They found a few 
emission lines (a broad O vii complex and C vi Lya line, plus 
a narrow Oviii Lya line), and several absorption features that 
they attribute to a multi-component WA. They fitted three xabs 
components to the combined RGS spectra obtaining ionization 
parameters fully in agreement with ours (log^ - 0.89;^q[j, 
2.14^qJ2, and 3.26^q27, respectively). However, there are inter- 
esting differences between the Smith et al. (I2007I I analysis and 
ours. First, we find a correlation between the ionization degree 
and the outflow velocity, i.e the more ionized is the component, 
the faster it flows. Smith et al. (2007) report the opposite, with 
their low ionization component having a faster outflow velocity. 
Further, their reported column densities for two of the three com- 
ponents are almost one order of magnitude higher than ours. It is 
difficult to ascertain whether these values correspond to intrinsic 
variations in the properties of the WA with respect to our obser- 
vation. Some probable explanations of these discrepancies are 
discussed in Detmers et al. (1201 Oi l. Given the fact that the ion- 
izing luminosity and the ionization parameter in the Smith et al. 
(2007) analysis are comparable to ours, we speculate that their 
higher column densities are due to their extremely low velocity 
broadening, which they find close to 0, probably because of the 
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Fig. 5. Qiandra-LETGS spectrum of Mrk 509. The solid line represents the best-fit model. The most prominent features are labeled. 



low signal-to-noise ratio of their data. In fact, their second WA 
component, for which they report a nns velocity of TO^^q km s"' , 
is fully consistent with ours. 

To further test this possibility, already suggested in Detmers 
et al. ( 1201 Oi l, starting from our best-fit values using a power- 
law as the continuum model (for consistency with the contin- 



uum model of Smith et al. 120071 ) we have set to zero the velocity 
broadening of our components 1 and 3 and re-fitted the data. 
This renders Nh - 2.0;^j 2 x 10^' cm"^ for the high-ionization 
component, which is in agreement with the value of Smith et 
al. (2007) at the Icr confidence level. The column density of the 
low-ionization component, however, does not change substan- 
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tially from its original value, remaining significantly different 
from of Smith et al. (I2007I I. Repeating this exercise with a spline 
model for the continuum we find a column density for the lowly 
ionized component which is consistent with that for the power- 
law model. 

Mrk 509 was observed in April 2001 for ~59 ks with 
Chandra-HETGS simultaneously with HST-STIS. Yaqoob et al. 
(12003 ) fitted the observed absorption lines using a single pho- 



;+0.39 
-0.45 



toionization model with logf - 1.76^g {^ and Nh - 2.06 
10^' cm"^. This best-fit has parameters between our components 
1 and 2. However, it has a larger column density (see Table |7]). 
Interestingly, they found this component to be outflowing at 
~200 km s"', in good agreement with the outflow velocity of our 
component 2. We note that this velocity was not a free parame- 
ter in their XSTAR modelling, but was introduced as a nominal 
velocity offset. The chosen value was based on the velocity pro- 
file analysis of the most relevant absorption lines present in the 
HEGTS spectrum. There are hints of a more complex kinematic 
structure since some high-ionization lines, such as Mgxi and 
Nex, show larger blueshift values (e.g. ~ -620 km s"' forMgxi) 
but the uncertainties make it difficult to ascertain whether there 
is an actual higher velocity component. Yaqoob et al. j2003l l 
also applied a curve of growth analysis to deduce a velocity 
width of 100 km s"', quite similar to what we find and fully 
in agreement with the recent analysis of the data gathered dur- 
ing our 2009 multi- wavelength campaign (Kaastra et al. 2011b) . 
hereafter Paper II. The velocity widths found in X-rays are, in 
general, broader than those found in the UV (Kriss et al. 120001 
Kraemer et al. 120031 ) which may provide some hints on the na- 
ture of the X-ray /UV absorbing gas. Note, however, that the UV 
absorption is characterized by more components at velocities too 
close to be resolved in the X-rays. This may explain the widths 
measured in the Chandra spectra if there are multiple X-ray ab- 
sorption systems corresponding to each of those seen in the UV 
(see Sect.l44b. 

Another set of three XMM-Newton observations were taken 
in 2005 and 2006. Cappi et al. (120091) used the pn spectra in 
these datasets and obtained evidence for a highly ionized and 
mildly relativistic outflow, while Detmers et al. (2010) used the 
RGS spectra to further analyze the WA in Mrk 509. They con- 
firmed the presence of three ionization components outflowing 
at distinct velocities, with log^ = 1.0 + 0.2, log^ = 1.95 + 0.02, 
and log ^ = 3.20 + 0.08, respectively, by stacking the three RGS 
spectra. The most remarkable difference with our results is that 
the outflow velocities measured by Detmers et al. (120101 ) are not 
correlated with the ionization parameter, finding that the inter- 
mediate ^ component is roughly at the systemic velocity of the 
source, while the low- and high-ionization components are out- 
flowing at -120 + 60 and -290 + 40 km s"', respectively. An 
analysis of the individual RGS spectra yielded similar results 
albeit with larger error bars. Interestingly, they found a hint of 
changes in the WA between the 2005 and 2006 observations, 
taken ~6 months apart. Unfortunately the statistics prevented 
them from constraining the location of the different WA com- 
ponents, except for the high-ionization one which was found to 
lie within 4 pc from the ionizing source. 

As part of our 2009 multi-wavelength campaign, Mrk 509 
was observed again by XMM-Newton for a total of 600 ks, split 
in 10 observations of 60 ks each taken 4 days apart in October- 
November 2009. The detailed analysis of the stacked RGS spec- 
trum is reported in an earlier paper in this series (Detmers et al. 
120111 hereafter Paper III). Although the sensitivity of the stacked 
600 ks RGS spectrum is much higher than our LETGS observa- 
tions, the results between the WA analyses of both observations 



are fully consistent within the error bars. The most relevant result 
is that this deep XMM-Newton observation finds no evidence 
that the ionized outflows consist of a continuous distribution of 
ionization states, but rather they are made of discrete compo- 
nents. Indeed, Paper III finds 6 ionization components (two for 
the slow outflow, consistent with the systemic velocity of Mrk 
509, and 4 for the fast outflow). If we compare our best-fit pa- 
rameters reported in Table Q with those of Table 5 in Paper III 
we clearly identify their component Bl with our component 1, 
whereas their component E2 corresponds with our component 3. 
On the other hand, our component 2 seems to be a blend of their 
components C2 and D2, which are unresolved in this Chandra 
dataset. 

Note, however, that the outflow velocities measured in this 
LETGS spectrum are consistently larger than those reported for 
RGS in Paper III for both xabs and slab fits. An analysis of 
the wavelength scale was carried out in Paper II, where a best-fit 
observation/model wavelength shift of -1.8 + 1.8 mA was re- 
ported combining all the indicators (see Table 1 1 in Paper II). In 
the same table, the difference in the wavelength scale between 
the RGS observation and the LETGS observation is quoted as 
-7.3 + 3.8 mA. Hence, there is a difference of ~5.5 mA be- 
tween the adopted wavelength in RGS and that of LETGS, the 
latter being larger (redshifted) with respect to the former At a 
typical wavelength of 20 A, 5.5 mA corresponds to a velocity 
offset of 83 km s '. Therefore this offset should in principle be 
substracted to the outflow velocities quoted in Tables |6] and Q 
if one wants to compare with the equivalent models reported in 
Paper III. However, there is, within the systematic uncertainties, 
no discrepancy between RGS and LETGS wavelength scales. 

Interestingly, the average flux state of Mrk 509 during the 
XMM-Newton observation was ~60% higher than during the 
Chandra observation in the soft X-ray waveband (20 - 35 A), 
and ~30% higher in the hard (7 - 10 A) band. Since the maximal 
variation happens at low energies we expect very little effect on 
most of the ionization states observed. The variations at harder 
energies may have an effect but the WA parameters remain sim- 
ilar in both observations. It is then difficult to state whether this 
is due to an intrinsic lack of variation of the gas (i.e. the gas did 
not have time to recombine after the flux decrease) or that there 
is an actual change in the parameters but it is statistically unde- 
tectable. Indeed, errors for some of the key parameters such as 
the column densities, which would allow us to put some con- 
straints on the location of the WA, can be as high as ~50%. 

4.2. Structure and location of the warm absorber 

In this section we discuss the nature and location of the ionized 
absorbers of Mrk 509. The low-ionization phase (component 1 
in Table |7]) appears to be slightly redshifted and is consistent 
with being at the systemic velocity at the 2cr level, while com- 
ponents 2 and 3 are outflowing at several hundreds of km s '. 
This suggests the possibility that component 1 is part of a dif- 
ferent structure than components 2 and 3. To test this we plot 
the pressure ionization parameter S as a function of the elec- 
tron temperature T. Components sharing the same H value in 
this thermal stability curve would likely be part of the same 
long-lived structure. The pressure ionization parameter is de- 
fined as H = L/4nr^cp = ^/6nckT. It was calculated using 
CLOUDY, which provided a grid of ionization parameters ^ and 
the corresponding electron temperatures for a thin layer of gas 
illuminated by an ionizing continuum. The resulting curve di- 
vides the r - S plane in two regions: below the curve heating 
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log S 

Fig. 6. Pressure ionization parameter as a function of electron 
temperature. The filled squares represent the warm absorber 
components in Mrk 509. 



dominates cooling, whereas above the curve cooling dominates 
heating (e.g. Reynolds & Fabian 1995 , Chakxavorty at al. 2008 ). 
Therefore, the parts of the curve where dT/dE < (where the 
curve turns back) are unstable against isobaric thermal perturba- 
tions (vertical displacements from the curve), while those with 
positive derivative dT/dE are stable. 

If we overplot the WA components on top of this stability 
curve (see Fig. |6]l we can see that the three of them fall in stable 
areas, which suggests that the absorbing gas belongs to a dis- 
crete scenario rather than being part of a continuous distribution 
of ionized material. Furthermore, the value of H for components 
2 and 3 are consistent within the error bars. This would mean 
that both components nominally share the same pressure ioniza- 
tion parameter, i.e. they are in pressure equilibrium. On the other 
hand, component 1 has a significantly different S value which, in 
addition to the fact that its features are not blueshifted, strongly 
suggests that it may not coexist as a part of the same outflow 
structure as components 2 and 3. This will be further discussed 
inSect.l4!4l 

The location of the WA is hard to estimate since the product 
nR^ is degenerate. The best way to put constraints on the distance 
R is therefore by using variability to determine the densities of 
the WA (e.g. Longinotti et al. 120101 Detmers et al. I2008I I. As 
mentioned in Sect. 14.11 there was a variation in flux of ~30% 
between the XMM-Newton 600 ks observation and the Chandra 
observation that would help to constrain the distance if changes 
in the WA were detectable. For that purpose we compared the 
measured column densities of the different ionization states of 
oxygen, as these ions were the most accurately measured in both 
datasets. At a first glance the ionic column densities measured 
in both cases are consistent within the error bars (see Table |6] 
and Table 4 in Paper 111). We statistically tested this by apply- 
ing a two-dimensional Kolmogorov-Smirnov test on both ionic 
column density distributions of oxygen in the same manner as 
in Ebrero et al. (,2010) . The test was performed for components 

1 and 2 in model 2 of Table |6] and its equivalent in Paper 111. 
The significance of the variation in component 1 between both 
datasets was at the ~90% confidence level, while for component 

2 it was well below the Icr level. Therefore, there are no grounds 
to estimate the distance of the WA absorber based on variability. 



Alternatively, assuming that the outflows are likely to be fil- 
amentary or clumpy, a volume filling factor defined in the same 
manner as in Blustin et al. (2005) can be characterized so that 
Cv < 1 . Therefore, if we integrate the density along our line of 
sight, we get A^h - CynR, where n is the density at the base 
of the absorber. Substituting n in the definition of the ionization 
parameter ^ we obtain 



^< 



LCy 



(1) 



We can estimate from our spectral fits the values of ^, A^h, and 
the ionizing 1 - 1000 Ryd luminosity L, whereas the volume 
filling factor can be estimated following dynamical arguments, 
assuming that the momentum of the outflow must be close to the 
momentum of the absorbed radiation plus the momentum of the 
scattered radiation (Blustin et al. l2005l l: 



C^ 



(f abs + f sc)^ 

1.23mpLv2Q' 



(2) 



where P^bs and P^c are the momenta of the absorbed and scat- 
tered radiations, respectively, nip is the mass of the proton, and 
Q. is the solid angle subtended by the outflow. We can estimate 
both momenta using formulas 20 and 21 in Blustin et al. ( I2005I I. 
and we take Q. = n sr, inferred from the type-l/type-2 ratio in 
nearby Seyfert galaxies of Maiolino & Rieke ( 119951 ). Using an 
ionizing luminosity of L - 3.2 x 10**^ erg s"' derived from the 
average SED of Mrk 509, and WA best fit values reported in 
Table |7] we obtain volume filling factors Cy of ~0.01, ~0.03 
and ~0.04 for components 1, 2 and 3 in the WA, respectively. 
Using now Eq.[T]we find that the outflowing components 2 and 
3 are located within 300 and 50 pc, respectively, from the central 
ionizing source. Using the relation of Krolik & Kriss (2001) the 
inner edge of the torus is estimated to be at ~5.5 pc from the ion- 
izing source. With the rather mild limits derived here is difficult 
to constrain the location of the WA, but it is consistent with the 
possibility that the outflows are generated in the putative torus, 
although an inner location cannot be ruled out. For component 
1, which is probably not part of the outflow, we derive a much 
larger upper limit of R < 6 kpc, giving support to the idea that 
this component probably belongs to the ISM of the host galaxy. 
In terms of cosmic feedback, we can estimate the mass out- 
flow rate for components 2 and 3 using Eq. 19 in Blustin et al. 
(l2005T l: 



Mv ~ Pabs + ^s< 



(3) 



In this way we obtain M2 ~ 1 Mq yr"' and M3 ~ 0.4 Mq yr"'. 
The total kinetic energy carried by the outflow per unit time (the 
sum of the kinetic luminosity of the blueshifted components), 
which can be defined as Lke = 1 /'2MoutV^, is of the order of 
4x1 0"*" erg s ' which represents a tiny fraction of the bolometric 
luminosity of Mrk 509, approximately LKE/^boi ~ 0.001%. It is 
then very unlikely that these outflows can affect significantly the 
local environment of the host galaxy. Indeed, current AGN feed- 
back models that reproduce the M-cr relation predicts that ~5% 
of the AGN radiative output is required to be fed back as both ki- 
netic and internal motion energy in the form of outflows in order 
to significantly affect the immediate surroundings of the central 
engine (e.g. Silk & Rees 1998 , Wyithe & Loeb l2003l Di Matteo, 
Springel & Hernquist 2005). The velocity widths that we mea- 
sure are of the order of 100 km s"\ approximately the same 
order of magnitude as the outflow velocity. The internal motion 
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Table 8. Best-fit parameters of the tentative ultra-fast outflow. 



-14,500^ 



Wavelength (A) 

Fig. 7. Chandra-LETGS spectrum of Mrk 509 in the ~ 5 - 
9 Aregion. The solid line represents the best-fit model. The most 
prominent features are labeled. Note the shift of the troughs with 
respect to their wavelengths at the rest frame of the source (ver- 
tical dashed lines). 



energy of the gas is hence similar to its bulk kinetic energy but, 
nevertheless, it is still far too low to generate a significant impact 
even if we consider a less conservative energetic requirements 
to disrupt the local ISM of Mrk 509 (Hopkins & Elvis 2010). 
Unless a higher velocity and mass outflows are present (albeit 
undetected in the LETGS, but see Sect. I4.3l l, the shallow and 
relatively slow outflows detected in Mrk 509 are therefore not 
powerful enough to play an important role in cosmic feedback. 

4.3. Is there an ultra-fast outflow? 

The region of the spectrum at shorter wavelengths, between ~5 
and 9 A, shows a noisy continuum although some absorption 
troughs could be present. The most prominent ones were those 
at -4.95 A and ~6.1 A. A tentative identification of these lines 
suggested that they could be S xv and Si xiv, respectively. If this 
was the case that would mean that the gas responsible for these 
features is outflowing at ~ 14, 000 km s"' . We tested this possi- 
bility by fitting a slab model in which the outflow velocity and 
the ionic column densities of S xv. Si xiv, Al xm, and Mg xii (the 
most prominent lines expected in this range of the spectrum) 
were the free parameters, while the widths of the lines were 
kept frozen to 100 km s '. The addition of this component im- 
proved the overall fit of the spectrum only at the 90% confidence 
level and yielded an outflow velocity of -14, 500;^4jy km s"' (see 
Table[8]l. The spectrum along with the best-fit model is shown in 
Fig.E] 

Ponti et al. ( I2009I I also reported the presence of a possible 
ultra-fast outflow in Mrk 509 based on the analysis of highly 
ionized Fe absorption lines present in the XMM-Newton EPIC 
spectra of 2005 and 2006. They obtained an outflow velocity of 

14,500 + 3600 km s"', which would be in agreement with 

our marginal detection. They predict, however, that the outflow- 
ing gas should be highly ionized (log ^ = 5.15), whereas the con- 
centration of the ions in our detection peaks at values of log ^ in 
the range 2.7-3. 1 and would be hardly detectable at much higher 
ionization states. The Qiandra-HETGS spectrum of Mrk 509 re- 
ported in Yaqoob et al. (i2003i) does not show significant absorp- 



lon 



log Mo 



411L. 



logr 



Mgxii 
Al XIII 
Si XIV 
Sxv 



16.57±0.34 

<16.7 
16.71±0.39 
17.65±1.20 



2.75 
2.90 
3.05 
2.80 



Notes. '"' Outflow velocity, in units of km s"'; 

**' Ionic column density, in units of cm"^; 

*^' Ionization parameter at which the ion fraction peaks for each ion, in 

units of erg s"' cm. 



tion troughs at these wavelengths and therefore we are unable to 
confirm the presence of such an ultra-fast outflow. 

4.4. Connection with the UV 

Prior to this campaign, Mrk 509 was observed in the UV do- 
main by FUSE (Kriss et al. 12000^ in 1999, and by HST-STIS 
(Kraemer et al. 12003! ) in 2001, the latter simultaneously with 
Qiandra-HETGS (Yaqoob et al. 2003 ). Their spectra showed the 
presence of seven and eight different kinematic components, re- 
spectively, with approximately half of them blueshifted and the 
other half redshifted with respect to the systemic velocity of the 
source. 

As a part of our multiwavelength campaign, Mrk 509 was 
observed with HST-COS simultaneously with Chandra-LETGS. 
The data reduction processes and discussion of the results of 
this observation are presented in a companion paper (Kriss et 
al. 2011). 

Similarly to previous observations, HST-COS observations 
find that the UV absorption lines in Mrk 509 present 13 distinct 
kinematic components, ranging from -408 to +222 km s"' (the 
minus sign denotes an outflow). For consistency, we list the UV 
components in this paper using the same numbering as in Kriss 
et al. ( 1201 11 1, component 1 being the most blueshifted and com- 
ponent 7 the most redshifted. Note that Kriss et al. (120 Hi ) num- 
bered the UV components using kinematic criteria, while for the 
X-ray components we have adopted a numbering criterion based 
on the ionisation state of the gas. In Table |9] we list the main 
properties of the intrinsic absorption features in the HST-COS 
spectra, such as the velocity shift with respect to the systemic 
velocity and the ionic column densities measured for the C iv, 
N v and O vi ions, as well as the possible LETGS X-ray coun- 
terparts. Note that the HST-COS observation was taken in the 
1160 - 1750 A spectral range. Since the Ovi feature falls out- 
side of this range, the data relevant to this ion are those of the 
previous FUSE observation of Mrk 509 (Kriss et al. .2000j . 

4.4.1. Kinematics 

Attending to their kinematic properties, one can see that the gas 
responsible for the UV absorption and the X-ray absorbing gas 
share some of the components, suggesting that they might be co- 
located. For instance, the UV components 1 and lb, outflowing 
at -408 and -361 km s"', respectively, may correspond with the 
X-ray component 3. We note, however, that the resolution in the 
X-ray domain is not comparable with that of the UV, so a one- 
to-one component identification is not always possible. Indeed, 
some of the X-ray kinematic components may be a blend of ve- 
locities and ionization states (see Paper III). Within the error 
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Table 9. Properties of the HST-COS intrinsic features. 



UV Component 


l-out 


A'civ" 


A'nv" 


^ovi"'^ 


LETGS" 


1 


-408 + 5 


31.2 ±1.5 


107.0 ± 9.5 


215.0 ±47.2 


3 


lb 


-361 ±13 


16.3 ±6.1 


17.5 ±1.6 


154.9 ± 39.0 


3? 


2 


-321+5 


136.3 ±41.7 


149.0 ± 10.7 


566.2 ±118.0 


3? 


2b 


-291 ±6 


128.9 ±41.4 


130.6 ±11.3 


1248.1 ±395.5 


2? 


3 


-244 + 5 


47.7 ±7.2 


88.6 ±8.1 


675.2 ±157.1 


2 


4a 


-59 ±5 


66.3 ±2.1 


93.6 ±6.8 


804.5 ± 387.2 




4 


-19±5 


250.0 ±11.0 


323.6 ±12.2 


8797.0 ±4120.3 




5 


37 ±5 


36.2 ± 16.9 


38.3 ±11.1 


683.5 ± 429.3 


1 


6a 


79 ± 12 


5.6 ±2.3 


20.6 ± 6.4 


518.8 ±398.5 


1 


6 


121 ±5 


12.3 ±6.4 


56.1 ±4.0 


3436.1 ±6985.0 


1? 


6b 


147 ±8 


17.6 ±5.3 


5.5 ±1.8 


104.5 ± 234.6 




7a 


184 ±6 


3.6±1.1 


21.0 ±2.9 


660.2 ± 176.7 




7 


222 + 6 


6.4 ±0.8 


23.5 ± 2.7 


667.4 ±1105.3 





Notes. *"' Outflow velocity, in units of km s"' ; 
**' Ionic column density, in units of 10'" cm"^; 
''■' Based on FUSE observations; 
** LETGS WA component counterpart. 



bars. X-ray component 3 may also be partly identified with UV 
component 2. On the other hand, there is a clear correspondence 
between UV component 3, outflowing at -244 km s"', and pos- 
sibly UV component 2b at -291 km s"', with the intermediate 
X-ray component 2. The UV components 5, 6 and 7 are inflow- 
ing, their radial velocity redshifted with respect to the source 
rest frame, similarly as the X-ray component 1 . Both UV com- 
ponents 5 and 6a, and X-ray component 1 are flowing close to 
the systemic velocity of Mrk 509, which may indicate that they 
move mainly transverse to our line of sight. 

Since the inflow velocities of the UV components 6, 6b, 7a, 
and 7 are of the order of hundreds of km s ' , and the redshifted 
X-ray component is not in pressure equilibrium with the oth- 
ers, strongly suggests that these components are not part of the 
outflow and may have an extra-nuclear origin. In fact, Kriss et 
al. (12011 ) speculate with the possibility that we are detecting 
high- velocity clouds in the interstellar medium of Mrk 509 host 
galaxy, similar to the Complex C in our own Galaxy (Thom et 
al. i2008i) . 

4.4.2. Column densities 

The total ionic column densities of C iv and N v measured in 
the HST-COS spectra summed over all the kinematic compo- 
nents are of the order of log A^ion (cm^^) ~ 15, so they are not 
likely to produce substantial X-ray absorption. Indeed, our X-ray 
analysis yields an upper limit for logA^NV < 15.4 (see Table |6]l 
which seems consistent with the ionic column density measured 
by HST-COS but, however, cannot be used to probe if the UV/X- 
ray absorbers have a common origin. On the other hand, the 
O VI measurement based on FUSE observations in 1 999 provided 
higher column density which therefore could be also detectable 
in X-rays. 

Making direct comparisons between the UV observations of 
O VI and the X-ray, even among the individual O vi components 
seen with FUSE is difficult because of the high level of satura- 
tion in most of these features, as pointed out by Kraemer et al. 
(I2003I I. The quoted column densities from Kriss et al. (I2000I I for 
the strongest components should be treated as lower limits. The 
deepest O vi absorption is seen in UV component 4, which is 
at rest at the systemic velocity of the source (see Table |9]l. No 
X-ray component is found at this velocity in our LETGS spec- 



trum. Other strong (and saturated) O vi absorption components 
are seen in UV components 3 and 6, which likely correspond to 
X-ray components 2 and 1, respectively. The latter is redshifted 
and is probably not part of the outflow, but the former shares 
the same kinematic properties as X-ray component 2 so they 
are likely co-located. We note that the FUSE and LETGS ob- 
servations were taken almost 10 years apart, and therefore varia- 
tions in the O vi column densities, although unlikely, are possible 
(Ki-issetallMD- 

This poses the problem of disentangling whether the gas re- 
sponsible for the UV absorption and the X-ray absorption are 
the same. Kraemer et al. (2003) claimed that none of the UV 
absorbers detected in the HST-STIS spectrum of Mrk 509 could 
be identified with the HETGS X-ray absorbing gas described in 
Yaqoob et al. (I2003I I. as the sum of all the ionic column den- 
sities of O VII spread over the different UV kinematic compo- 
nents only accounted for the 10% of the column density for the 
same ion measured in the HETGS spectrum. The total value re- 
ported in Kraemer et al. ( I2003I I is logA^ovii - 16.51, which is 
indeed only 25% of the A^ovii column density we measure in 
our LETGS spectrum assuming that this feature originates in a 
single kinematic component (model 1 in Table |6]l as in Yaqoob 
et al. (120031 1. When we consider a multi- velocity wind (model 
2 in Table |6ll then the O vii column densities measured in ei- 
ther kinematic components are in agreement with the HST-STIS 
result within the error bars. Likewise, if we follow the same ar- 
gument considering the total O vi column density measured by 
FUSE, which sums up to log A^ovi = 16.1 and should be consid- 
ered as a lower limit, we find that this result is fully in agreement 
with ours when we assume a single kinematic component. More 
interestingly, the O vi column density in the UV component 3 
(logA^ovi = 15.4, also a lower limit) is consistent within the error 
bars with the measured column density in LETGS for the corre- 
sponding X-ray kinematic component (log A^ovi - 15.7 + 0.4). 

We see that the absorbing gas scenario in Mrk 509 is quite 
complex. The results obtained in this campaign show that the UV 
and X-ray absorbers have consistent column densities. Kriss et 
al. (120111) shows that the UV absorbers, although characterized 
in discrete kinematic components, are likely to be smoothly dis- 
tributed over a range in velocities, in particular for highly ionized 
species. If the UV- and X-ray-absorbing gases are indeed co- 
located, at least for some of the detected components, this could 
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be consistent with a clumpy scenario in which high-density low- 
ionization UV-absorbing clouds are embedded in a low-density 
high-ionization X-ray gas, as proposed in Kriss et al. (2000) and 
(1201 11 1. We defer a more exhaustive analysis on the relation be- 
tween the X-ray and UV absorbing gas including all the obser- 
vations available from this campaign (Qiandra-LETGS, XMM- 
Newton RGS, and HST-COS) to a forthcoming paper (Ebrero et 
al., in preparation). 



5. Conclusions 

We have presented here the results from the 180 ks Chandra- 
LETGS spectrum of Mrk 509 as part of an extensive multi- 
wavelength campaign. We found that the ionized absorber in- 
trinsic to the source presents three distinct degrees of ioniza- 
tion (log^ =1.06, 2.26, and 3.15, respectively) and is some- 
what shallow (A^H - 1-9, 5.4, and 5.9 times 10^° cm"^, re- 
spectively). The low-ionization component is slightly redshifted 
(Av = 73 km s"'), whereas the others are blueshifted with re- 
spect to the systemic velocity of the source (Av = -196 and 
-455 km s"', respectively, see Table [T]), and therefore they are 
outflowing. Their outflow velocities are correlated with the ion- 
ization parameter so that the higher the ionization parameter the 
faster is the outflow. The redshifted component is not in pres- 
sure equilibrium with the outflowing components, thus indicat- 
ing that it may not form part of the same structure. We put mild 
constraints on the location of the outflowing gas, which would 
be consistent with a torus wind origin, although an inner origin 
cannot be ruled out. The energetics of the outflows indicate that 
they have a negligible effect in the interstellar medium of the 
host galaxy. 

Mrk 509 was also observed in the UV domain with HST- 
COS simultaneously with Chandra, showing 13 kinematic com- 
ponents in the intrinsic UV absorption features. We compared 
the X-ray and UV absorbers in order to ascertain whether they 
have a common origin. At least three of the UV components 
could correspond with the X-ray components detected in the 
LETGS spectrum based only on the kinematics of the gas, which 
points to a possible co-location of both absorbing gases. The col- 
umn density of O vi measured in one of the X-ray components is 
in agreement with that of a previous UV observation with FUSE 
for their corresponding kinematic component. Similarly, the up- 
per limits on the column densities of C iv and N v found in X- 
rays are also consistent with the values obtained by HST-COS. 
Therefore, the similar kinematics and column densities indicate 
that the gases responsible for the X-ray and UV absorption share 
a common origin. 
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